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Abstract 

We continue our study of Horava-Lifshitz type theories using the methods of the spectral 
geometry. In this work we construct the infrared action of gravity and matter coupled to gravity 
in the most general way respecting the foliation preserving diffeomorphisms. This is done with 
the help of the spectral action principle based on some generalized Dirac operator. The gravity 
part reproduces the infrared limit of the Horava-Lifshitz gravity, while the matter part gives the 
generalization of the earlier suggested models. Due to the fact that the same Dirac operator 
is used in the construction of both sectors, the parameters of the gravity and matter parts are 
related. We expect that this potentially could naturally exclude fine tunings needed to get some 
desired properties as well as open new possibilities for the experimental tests of the model. 


1 Introduction 

In view of the ongoing search for the theory of Quantum Gravity and taking into account that 
neither String Theory nor Loop Quantum Gravity could not yet be considered as a completely 
satisfying solution, any attempt in this direction should be welcomed. One of the recent proposals 
is the so-called Horava-Lifshitz (HL) gravity [1]. This theory is based on the foliation-preserving 
diffeomorphisms (FPDiff) instead of the full diffeomorphisms of space-time. While the proposal 
has several interesting features expected from a consistent theory of quantum gravity, such as the 
power-counting renormalizability and ultraviolet spectral dimension of two m, it suffers from sev¬ 
eral problems (for the detailed account on the HL gravity and some of its applications, see [&H7]). 
One of the major drawbacks of the proposal, from our point of view, is a very large number of free 
parameters, which is of order of 100 [§]. To reduce this number one has to impose very restric¬ 
tive conditions like the detailed balance condition [T], which are neither natural nor satisfactory 
experimentally (see, though, m- 


* danielvieiral@gmail.com 

' arthurmamiya@gmail.com 

* apinzul@unb. br 


1 



Another, somewhat related, point is the coupling of the gravity based on FPDiff to matter. 
There are some arguments [10] that in contrary to the gravity part the matter sector should respect 
the symmetry under full diffeomorphisms (or it should be extremely fine-tuned). Clearly this is very 
non-natural and one would like to have a more symmetric formulation of the theory where both, 
gravity and matter sectors, are based on FPDiff or one should give more fundamental reasons why 
matter should respect much larger symmetry. 

In this paper, we are advocating the first possibility, i.e. that both sectors have FPDiffs as the 
fundamental symmetry. This is achieved via the spectral action principle mi- The main idea is 
that both sectors are controlled by the same object - some physically relevant Dirac operator. Due 
to this, both sectors are not independent and the parameters of the gravitational part are related to 
the parameters of the matter sector. Earlier this strategy was successfully applied to the Standard 
Model H2j. We are applying this approach to construct the most general infrared (IR) action of 
the HL gravity coupled to fermionic matter. The action is based on some natural generalization of 
the standard Dirac operator. In our earlier works mm we used the same type of Dirac operator to 
analyze the spectral dimension of the HL space-time (which was confirmed to be 2 for an arbitrary 
curved space-time) and to study the geodesic motion of a test particle in the HL gravity. The main 
result of the analysis of the geodesic motion was that in the case of the non-minimal coupling it 
differs from the motion naively calculated from the underlying Riemannian metric. Of course this 
result was not a big surprise, but what is important that in our approach the deviation from the 
Riemannian geodesics is controlled by the same Dirac operator that is used in the construction of the 
HL modification of the Einstein-Hilbert action. Potentially this might lead to some cancellations 
of the experimentally problematic effects. This urgently requires further studies. In the current 
work, we take this proposal to the next level and study the most general IR action in both, gravity 
and matter, sectors. Here IR means energies much lower then some typical quantum gravitation 
scale where the higher order terms in the HL gravity become relevant. Of course, this includes the 
range accessible to the present experiments and well beyond. The main results of the paper are the 
following: 

i) On the gravity side: The most general IR action of the HL type is derived as the spectral 
action for some generalized Dirac operator. The non-trivial result is that to obtain the most general 
form of the gravity action, one has to add P-violating terms to Dirac operator. 

ii) On the matter side: Using the same Dirac operator as in (i) and the matter part of the 
spectral action, we obtain the most general Lorentz violating fermionic action based on FPDiff and 
consistent with the Standard Model Extension (SME) of [14] , In our approach, the free parameters 
of SME are expressed in terms of the parameters of the Dirac operator. A very important result 
is that these parameters also depend on the geometry of the space-time, which can make Lorentz 
violating effects very small in a weak gravitational held. 

iii) The combination of (i) and (ii) naturally leads to the non-trivial relations between the 
parameters of the gravitational and matter actions. 

The plan of the paper as follows. In the section [2] we construct the geometric or gravitational 
part of the spectral action and compare it to the IR limit of the HL gravity. The section [3] deals 
with the fermionic or matter part of the spectral action. The following section discusses the possible 
phenomenological consequences of our approach. We collect all the notations and conventions, as 
well as the major part of the calculations in several appendices. 
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2 Spectral action: gravity part 


In this section, after the brief review of the spectral action principle, we derive the gravity action 
based on some generalization of the standard Dirac operator. We show that this action is exactly 
the IR limit of the general HL gravity. 

The spectral action principle m deeply roots in non-commutative or spectral geometry m- 
The main object in this approach to geometry is the so-called spectral triple. This is strongly 
motivated by the reformulation of Riemannian (i.e. Euclidean) geometry of a compact manifold 
M in terms of the commutative spectral triple: the commutative algebra C°°(M ) acting on the 
Hilbert space of square integrable spinors ^(S) and the standard Dirac operator. In this case it was 
possible to prove the so-called reconstruction theorem |16| that proves the complete equivalence 
between the standard and spectral approaches. The non-commutative formulation, being purely 
algebraic, opens up the possibility for various generalizations of the commutative geometry. This 
is achieved by generalizing in some consistent way one or all of the elements of the spectral triple, 
see m for some examples of the generalized geometries. Physical interest in the methods of non- 
commutative geometry exploded after the paper m but even before this there had been a constant 
interest in physical applications. For the purposes of our work the paper EU is most important. In 
this paper the concept of the spectral action principle was introduced and later it was successfully 
applied to produce a new approach to the Standard Model, see e.g. |12j . 

In words, the spectral action principle can be formulated as follows: the spectral triple contains 
not only the complete geometrical but also physical information. The main problem is to choose the 
right spectral triple. If the choice is made, one postulates that the dynamics of a physical system is 
governed by the following action 

S = Tlf + {J'lp-.O'l/j) = Sg eom + Smatt , (1) 

where O is some (generalized) Dirac operator, / is some cut-off function, A is some characteristic 
scale and J is a real structure (i.e. we should consider a real spectral triple |18j). While the S geom 
part contains the description of gravity and gauge fields, the S ma u term describes the natural and, 
in some sense, minimal interaction between the fermionic matter and the geometry (including gauge 
fields). As we have already mentioned, for the right choice of the real spectral triple the Standard 
Model perfectly fits into this scheme m- A very important point about m is that the same 
object, the generalized Dirac operator ID) controls both parts of ©■ As we briefly discussed in the 
introduction, this is the source of the possible relations between the parameters of two seemingly 
unrelated sectors. In this section we deal with the S geom part and discuss S ma tt i n the next section. 

To evaluate the trace in Tr / one has to use the heat kernel techniques. Here we just 

present the needed results, for details, see e.g. m- The analogous calculation for the standard 
Dirac operator was first done in (201 [21] , Because these methods (as well as the non-commutative 
approach to geometry) are well defined for Euclidean case, we will be working with Euclidean version 
of gravity. (This corresponds to the choice e = 1 in the calculations from the appendices.) Also we 
will work with a manifold without boundary. 

So, let P be some (generalized) elliptic operator of the order rn represented on a Hilbert space 
H of square-integrable sections of some bundle over d-dimensional manifold M. The heat kernel of 
P is defined by 

K(t,P) =Tr n e~ tp . 
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This heat kernel has the well known asymptotics for small t 

AT(i,P)~^r^Ta n (P) , 

n> 0 

where a n (P ) are completely defined by some local densities, a n (x,P), known as Seeley-DeWitt 
coefficients: 

a n (P) = [ a n (x,P)yfgd d x . (2) 

Jm 

There are several techniques how to calculate these coefficients one of the most effective being due 
to Gilkey [22] . Below we will give the explicit expressions for the first two coefficients for the case 
relevant for our purposes. The relation between the heat kernel and the trace of some function of P 
can be established through the operator analog of the Mellin transformation. We write this relation 
for the case of our interest, i.e. when d = 4 and m = 2. The fact that we work with a manifold 
without boundary can be used to show that a n (P) = 0 for all odd n [19]. Then the trace of /(P) 
is given by 

Tr/(P) = ^/ 2fc a 2fc (P) , (3) 

k> 0 

where f 2 k are 

oo oo 

fo = j f (u)udu, f 2 = j f (u)du, f 2 {k+ 2 ) = (-l) n / (n) (0), k > 0 . 

0 0 

In m it was shown that choosing P = D 1 2 where D is the standard Dirac operator on the 
twisted spinor bundle over M, the first three terms in ([3]) reproduce General relativity with the 
cosmological constant coupled to some gauge field (depending on the twisting). 

The idea of this section is to apply m to P = D 2 , where D is some physically motivated 
deformation of D. As it was argued in mm, the most natural Dirac operator that respects the 
fundamental foliation structure of the HL geometry should be written as the most general differential 
operator represented on 'H = L 2 (S) such that it is of the first order in time derivative and up to 
the third order in space derivatives. Each term should have well-defined transformation properties 
under FPDiff. To write such an operator explicitly, the (3+1) decomposition of the standard Dirac 
operator was worked out in details in m 

D = = 7 °D n +( 3 )p-I 7 °A- + . (4 ) 

Here D n is roughly speaking the covariant derivative along the normal vector, ( 3 (D is Dirac operator 
on the 3d leaf of the ADM-like foliation of space-time, K - its extrinsic curvature and N is the 
lapse function. (For the details of the derivation and used notations see [13] and Appendix lAl of the 
present paper.) While the general high order generalized Dirac operator is very hard to studyj^] its 

1 For example, even the study of the heat kernel of the flat case of the third order generalized Dirac operator is 

already technically quite involved problem [23| . Nevertheless, in |3] it was shown that some conclusions based on the 

most general form of the generalized Dirac operator suitable for a HL-type theory can be made. In particular, it was 
shown that the UV spectral dimension will always be equal to 2. 
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IR limit, i.e. when only the lowest order derivatives are kept, can be analyzed. In [T3j j we argued 
that the most general form of the deformed Dirac operator in IR limit is 

D = 7 °D n + ci ( %+c 2 7 °/i + c 3 7 “a Q + c 4 it + c 5 7 ° 7 a a a , (5) 

where q, i = 1, 5 are some arbitrary parameters and we introduced the notation a a = Actually 
in m we used ([5]) with c 4 = C 5 = 0. While this does not affect any of the conclusions of that paper, 
as we will see, adding these terms makes the real difference for the present work. (|5]) is exactly 
the operator that we want to use as the fundamental object in (JT]). One can bring it to a more 
convenient form. One of the conclusions of m was that from the point of view of the geodesic 
motion of a test particle, the original metric on M, g )W , is substituted by some effective one with 

the re-scaled space part, g uu , see Eq. (140l) of the appendix IU1 with a = \. In the appendix iDl we 

^ f 1 

show that in terms of this more natural metric ([5|) can be re-written as (154[) 

D = 7^ V“ + ^c 2 + 0 K' y° + ^c 3 - y7+ c A K + c 5 j 0 fe/ l a IJi . (6) 


So it takes the form 

D = + F , where F = ^C 2 + 

To calculate the geometric part of the spectral action (JT]) we need to calculate D 2 . In the appendix 
m we show that one can derive the generalization of the Lichnerowicz formula for the generalized 
Dirac operator of the form (0 


7 ) K l'° + ( c 3 - ) f e t IJ ' a p. + C4 K + C57°7*e i /i a At . ( 7 ) 


D 2 = + E ? ( 8 ) 

where D and E are some connection and endomorphism respectively, see eqs. (1361) and (1371) of the 
appendix iBl evaluated for the re-scaled metric (1401) . The importance of the representation (JHJ) is 
that there is the general procedure of calculating the Seeley-DeWitt coefficients that enter in (J3|). 
For our purposes we will need only the first two (T9li22 j 

a 0 (x,D 2 ) = (47r)~ 2 Tr(l) , 

a 2 (x,0 2 ) = (4vr)- 2 TV^-|l + ^ , (9) 

where we already set that d = 4 and the trace is over the spinor indices, i.e. in our case Tr(l) = 4@ 
Here R is the 4-d curvature calculated for the re-scaled metric (1401) . I 11 the appendix [0 we show 
how it can be written in terms of (3+1) decomposition, see (1491) 

R = c\ ^R +K 2 - K I _ W K 111 ' + 2V m ( c\n x V x n^ - n^V A n A ) . (10) 


2 The dependence on Cl will appear only in a+a;,© 2 ) but the contribution of this coefficient will be suppressed by 
factor. 
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The trace of the endomorphism E is calculated in the appendix [B] (1391) . To apply this result to our 
Dirac operator d7j) we have to identify the coefficients a, b^ and c^ v from the appendix formula ()38[) 
with the parameters in (0: 


a = caK , b a = 


c 2 T ~ ) K 


C 3 


Cl 


LI 

e i a » 


; CQi — C5 G^CLfj, , Cij — 0 


where we passed from the space-time indices to the flat (tangent) ones. Then (1391) will become 

Tr E = — 12c\K 2 + 4c 2 a M a M + R . (11) 

Using COD and (fill) in 0 we can easily get the integrated Seeley-DeWitt coefficients 0 


«o(B 2 ) = ^2 J Vdd A x , 

a 2 (B 2 ) = —^ J ( c ? (3) ^+(l-36cl)K 2 - + 12c|a%) y/gd*x . (12) 

Finally, using (fl2l) in Q we arrive at the expression for the geometric part of the spectral action 

Sgeom 

Sgeom = j ^ C 1 (3) A+U ~ 36cf)AT 2 - + 12^0^ + \f9 < £x . 

This should be compared with (the Euclidean version of) the IR limit of the action for HL gravity 
with the cosmological constant A c |241l25| 

i\/f2 p , v 

SiRHL = ^f U { 3 ) R+^K 2 - K tJ ,uK> w + 770% + A c ) ^d 4 x . ( 13 ) 

2 ijl v ' 

It is obvious that upon the following identifications of the free parameters: 

M2p= 2&' ^T = e, A = (l-36c 2 ), r ?=12 c|, A c = (14) 

both actions are exactly the same. This is the main result of this section. Using the map between the 
parameters (1141) we can translate the experimental bounds on (£, A, 77 ) into the bounds on (ci, C4, C5), 
see below. It is important to note that without the P-violating terms in 0, which are proportional 
to C4 and C5, the resulting action would be a rather trivial deformation of GR. Also note that up to 
this order there are no bounds on ( 02 , 03 ). This should be compared to the same conclusion in [13]. 
This situation changes when we include matter. So we pass to the consideration of the matter part 
of the spectral action 0 , S ma tt- 


3 Spectral action: matter part 

The spectral action principle 0 tells us that the action for fermionic matter is given by 

Smatter = (Jlp = J d 4 X . 


(15) 
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Before we proceed with our analysis of (1151) some comments are in order. The matter in HL 
gravity has been considered before with the conclusion that to avoid problems (fine tuning, strong 
Lorentz violation, etc) one has to consider the minimal coupling of matter sector to the geometry, 
i.e. the same as in GR [10l[26j . We can immediately see that this is against the spirit of the 
spectral action principle - clearly the matter in (1151) is not minimally coupled in this sense. But it is 
minimally coupled from the point of view of the spectral action - the action (fl5l) is really minimal 
and the most natural action built of some physical Dirac operator. The fine tuning, as we said 
before, could be if not avoided completely but at least improved by the fact that the same Dirac 
operator controls both parts of the full action (see also the discussion below). 

Recalling the form of our generalized Dirac operator (j5]) it should be obvious that the action ()15[) 
will lead to Lorentz violation, that can in principle be detectable in experiments, which, in its turn, 
will put severe constrains on the parameters in ([5]). In |14| . a general model independent Lorentz 
violating extension of Standard Model (SME) was presented and in [57] . the coefficients appearing 
in | I14| are compared to precision experiments. Our goal here is to see how (fl5l) fits into this scheme. 
As we have stressed many times the same free parameters c* enter both the geometric (gravitational) 
and fermionic (matter) parts of the action. So if we take the spectral action principle seriously, then 
we can use the data from cosmology or gravitational precision measurements to constrain Standard 
Model deviations from Lorentz invariance, and vice-versa. Of course to do so one has to work with 
the complete SM (with all the gauge field content). We will use our (somewhat toy) model to 
demonstrate how this should work in principle. 

To proceed, let us write down the most general SME action in the fermionic sector [H1J 

S = J d 4 x^(e a ^T a V^ + ijjM'ip) , (16) 


where 

pa 

M 


r - C/H ,e°V 7 ft - VtV - e ^ - if 
m + inn 57 s + m^'f + V/tV + 


1 


av 


e fe V7 &c 


(17) 


Because (|16[) is model independent, all the parameters in (1171) are independent. If one tries to arrive 
at (1161) within a specific model, typically there will be some conditions on these parameters. E.g. 
this happens when one re-writes the noncommutative electrodynamics in SME form, see |28( . To 
see how this works in our case we have to re-write (1151) in the form (I16|) . Here, actually, we have 
two alternatives. 

i) Using (1521) . we can easily show that the Dirac operator ([5]) can be written in the form 


D = + ( Cl - l)e/ 7 ‘V; + (| + c 2 ) 7 0 A" + 



7*eAa M + c 4 K + . 


(18) 


ii) Alternatively, we can directly use the representation Q 

J5 = rv - + ^ C2 + 0 RT7 0 + ^c 3 - y^ 7*e+ c^K + . ( 19 ) 

Naively, it would seem that (fl8l) is the right choice because its first term represents the standard 
Dirac operator with respect to the original metric g^ u . But taking into account the conclusion of the 
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paper [13] that the physical geodesic motion of a test particle is defined as the usual geodesic motion 
but with respect to the re-scaled g^ u metric, the choice (1191) might seem more natural physically. 
In any case, this requires some further understanding and here we treat both cases. 

i) Comparing the first choice (|18[) to ([16I17D we easily identify the non-vanishing coefficients 


Cfii/ — (Cl l)^htiy 


rrif, = 



K 



m = c^K , Hq 0 = 2c§a a 


( 20 ) 


To get we used (l30l) . 

ii) Doing the same for the choice (fl9l) we immediately arrive at the following set of the Lorentz 
violating coefficients: 


mu = 



K 



, m = c 4 K , H 0a 




( 21 ) 


the rest of the coefficients being zero. 

So the difference between two cases is minimal and essentially reduces to the absence of the c„„ 
correction in the second case, which is predictable due to the metric re-scaling. 

There are several short comments that we would like to make before proceeding with the dis¬ 
cussion of the obtained results. 

Firstly, the matter part of the spectral action finally gives some non-trivial dependence on the 
coefficients C 2 and C 3 . As we have seen in the section [ 2 ] these coefficients completely drop out of 
the geometric action. This, of course, agrees with the earlier observation that the corrections to the 
physical geodesic motion due to the generalized nature of ([5]) depend only on c\. 

Secondly, it is important to notice that all the Lorentz violating coefficients depend on the 
background geometry. This is very important when comparing to experiments because the absence 
of the measurable corrections would not automatically mean that they are small in general but could 
be just because the gravity is weak in this particular experiment (see also below). As an example of 
such correction we can give the gravity dependent mass: even if the fermionic field has a vanishing 
bare mass m = 0 , the term 


Sjnass = j (fxy/g^CiKll) 

will generate an effective mass different from zero, dependent on the extrinsic curvature of the space. 

Thirdly, from (l20l) or (1211) it is clear that, of course, the natural way to parameterize the Lorentz 
violating sector is not by c* but by the deviations of these coefficients from the undeformed values 
{d} = (1, —1/2,1/2,0, 0) (this can be easily derived by comparing (J4|) and ([5])). 


4 Constraining the parameters 

Typically the bounds on the parameters of Horava-Lifshitz gravity are coming from the experimental 
tests where gravitational effects are significant. Usually, the experiments in principle capable of test¬ 
ing the theory do not have enough precision do differentiate from usual GR or are still unrealizable 
(see for example [2411251129] and references therein). If however the spectral action formalism is the 
correct one to deal with HL gravity, than the parameters of the gravitational action are related to 




the parameters in the fermionic action and those are constrained by much more feasible experiments. 


As described in the previous section, in the spectral action formalism the parameters of the 
effective IR theory for the bosonic and fermionic degrees of freedom are both obtained from the same 
object, the anisotropic Dirac operator. When we choose the Dirac operator j5J) we reproduce, to the 
lowest order, the IR limit of Horava-Lifshitz theory and a Lorentz violating fermionic action. In m, 
the general Lorentz violating extension of QED was presented (11611 . and in m the experimental 
bounds for the general parameters in (22) were discussed. These parameters are related to the 
parameters we introduced for the Lorentz violating Dirac operator (11811 or (11911 via (12011 or m 
respectively. The same parameters reproduce the Horava-Lifshitz action (to the leading order) via 
G3D- Thus, combining (1141) and (12011 . we gct(^| 


C/iv — 


(£ 2 


1 


= 



K 



m 


= -Vl^XK , H 0a = Jl 


( 22 ) 


Therefore, the experimental bounds on the particle side relate directly to the gravity side. 

We can rewrite (122(1 using parameters (/3, fx, a ) more convenient to compare with astrophysical 
data [25] 

c _ 1 x ol 

s 1 3 '' 1-/3 ’ ^ 1-/3 ’ 

so (m a is unchanged) 


(£-l) ; 


- 1 


hfxv i n^o — 


L 2(/3 — l) 2 


+ c 2 


IV 1 u I — 

K , m = -\ - -K , Ho a = . /--< 

’ 6 y 1 — /3 V 3 - 3 / 3 


(23) 

It is instructive to write the relations above as a table connecting the Lorentz violating param¬ 
eters as in m- the parameters of the anisotropic Dirac operator (fT8l) q or the IR Horava-Lifshitz 
coefficients (£, A, 77 ) and the ADM functions (lapse function, 3d metric and extrinsic curvature): 


LV parameter 

HL/anis-Dirac parameters 

ADM functions 


£ (or /3) 


m 0 

f (or /3), c 2 

K 

m a 

C 3 

d a 

m 

A (or /3 and /i) 

K 

Hq & 

r] (or /3 and a) 

ci a 


The spectral action formalism provides us a connection among the three columns. The first 
column exhibits Lorentz violating parameters for the action describing matter (fermionic) fields. 
Experiments testing the propagation and dynamics of fermions can provide bounds on these pa¬ 
rameters, even without gravitational effects. In [27] it is presented many such experiments and 
bounds. 

3 If we instead use the Dirac operator with the rescaled metric m we will get for (1231) below 

To avoid cluttering in the following analysis we consider only the first case, but all the arguments follow analogously. 


3-3 $ 
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The second column contais the parameters that measure deviations of HL gravity from standard 
GR, and are therefore constrained by gravitational experiments, i.e., experiments where the gravita¬ 
tional field is significant. For example, in 1 241125] it is presented how measurements of the radiation 
emitted by the accretion disks around black-holes or gravitational lensing can provide bounds on 
the HL parameters. In [29] Solar System tests (radar echo, perihelion precession of Mercury and 
light deviation by the Sun) are presented. 

The third column is related to the gravitational field and the foliation structure of spacetime and 
thus, in principle, can be obtained directly from theoretical calculations. For example, in appendix 
El we show how to write a a in the next-to-leading order in Horava-Lifshitz parameters for the case 
of a weak field for earth-based experiments and obtain the explicit form for the radial component 
of the Lorentz violating parameter m^ f)55[) 


m r 



G n M 

^2 



g 2 n m 2 v 

4 r 2 £ 


As an example of how to connect the three columns presented above, consider the experimental 
constrain on the HL parameter j3. In [241125) . it is argued that by analyzing the frequencies and 
redshifts of X-ray emitted by accretion disks of black-hole candidates]!! their emission lines and 
the impact parameter of the circular photon orbit (that affects gravitational lensing) one can in 
principle measure possible Horava-Lifshitz deviations from the standard Einstein gravity. However, 
these deviations from GR are too small to be conclusively measured (at least at the present time). 
For instance, a large value of (5 = 0.2 (that is, £=1.25) gives a deviation of the redshift measured in 
accretion disks by less than 0.1% when compared to the GR prediction. 

If however we take the spectral action formalism as the correct way to deal with effective actions, 
we can instead look at experiments in the matter sector where we have much better precision. 
Altschul [30] gives bounds on the electron Lorentz violating parameters cjj by analyzing synchrotron 
radiation and inverse Compton scattering from astrophysical sources, with I, J = X, Y, Z. Here the 
coordinate basis is chosen in a way that Z corresponds to the Earth north, i.e. parallel to Earth 
rotation axis, X corresponds to the direction from the Earth to the Sun at the vernal equinox, and 
Y follows the right-hand rule. E.g., the bound on cxx is given by 


- 5 x lO" 15 < cxx < 5 x 10" 15 . 


In the spectral action formalism cxx is related to f3 via (123[) 


(24) 


Cxx 


1 


hxx ■ 


These electrons are produced in different astrophysical sources, so for each one we would have 
to take into account hxx in order to calculate /3 . If we consider sources where the gravitational 
field is weak, approximating hxx ~ 1 we get 


|/3| < 2.5 x lO’ 15 , (25) 

which is a much stronger bound then the one coming from any of the gravitational experiments. 
Thus, while the still challenging gravitational experiments may not provide us with bounds on GR 
modifications, other kinds of high precision experiments may. 

4 GR modifications change the Innermost Stable Circular Orbit (ISCO), and therefore modify accretion disks. 
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5 Discussion and conclusion 


In this work we studied the application of the spectral action principle (H|) to the gravity plus matter 
theory based on the foliation preserving diffeomorphisms. Choosing © as the physically motivated 
generalized Dirac operator in the infrared regime, we were able to reproduce the most general IR 
action for HL gravity coupled in a non-trivial way to fermionic matter. 

One of the main conclusions that can be made based on our results is that, if taken seriously, 
the spectral action principle can serve as an effective tool in reducing the parameter space of the 
theory. E.g., in the IR regime, considered in this work, the parameter space is just 5-dimensional 
one, spanned by {cj}- 

Another very important (and potentially could be the most important) advantage of our ap¬ 
proach, intimately related to the previous one, is the possibility to use some high precision experi¬ 
ments from one sector (typically this would be the matter sector) to constrain the parameters of the 
other one. We demonstrated how this could be done in principle on the example of the parameter 
/3, which does not have very restrictive bounds coming from gravitational experiments, see (1241251) . 

We have not done any serious confrontation of our model with the existing experimental data. 
This is because our model is still a toy model that served well to our purpose - demonstration of the 
advantages and the potential of our approach. Of course, the most obvious and urging next step is 
to make our model less a toy, which will allow to look at the experimental bounds on parameters 
of the model. This will require adding the gauge sector to our construction. The spectral action 
principle provides a natural way for the inclusion of gauge fields - they will become a part of a 
generalized Dirac operator on some appropriately twisted spinor bundle. The example for the usual 
case is considered at the end of the appendix[B] For the case of a theory based on FPDiffs, one has to 
construct the twisting consistent with this symmetry. This problem is currently under consideration. 

As we stressed several times, our approach could help to avoid some fine tunings needed, as it 
was argued in [lOj , for the consistent non-minimal coupling of matter. To check whether this is the 
case, one has to do the analysis along the lines of m but now with the constrained parameter space. 
The optimistic expectation would be that some problems observed in |10] will be cured by the built 
in relations between the parameters. We are planning to address this problem in the nearest future. 

Much more ambitious and fundamental problem is the generalization of our approach to the full 
theory, i.e. not just to the IR limit. As we commented above this would require working with the 
most general deformation of Dirac operator by the terms up to the third order in space derivatives. 
This is a very challenging, both technically and conceptually, problem. Some initial steps were 
taken in [23], but a lot of work still should be done before one could try to derive the full HL theory 
from the spectral action principle. Taking into account the advantages of this approach, we still 
hope that this will be possible to accomplish. 
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A Notations and conventions 


Here we fix the notations and conventions used in the main text as well as establish some formulas 
used in the further calculations. Because this is very standard material, we essentially copy one of 
the appendices of [13]. Though we will need only the Euclidean signature, for the sake of generality 
we will obtain our results for both, Euclidean (e = 1) and Lorentzian (e = —1), cases. So the only 
difference from |13| is that now we explicitly write all formulas for both signatures and use the 
algebra of 7 -matrices that differs by a sign from m- 

Coordinate system. We adopt the coordinates compatible with the foliation structure: 

= ( t, x ) , (26) 

where t = const defines a leaf of the foliation E^, while x are the coordinates on E*. 

Metric. In the coordinates (126[) . the metric takes the usual ADM form m 

ds 2 = e(Ndt) 2 + h a p(dx a + N a dt)(dx 13 + N^dt) , (27) 

where N is the laps function and N a is the shift vector and h a p is the 3d metric on a leaf. Throughout 
the paper we are using the following system of indices: 

The Greek letters from the middle of the alphabet, /a, v, ... are used to denote the curved 
coordinates ()26l) and take values 0, 1, 2 and 3. 

The Greek letters from the beginning of the alphabet, a, (3, ... are used to denote the space 
part of the curved coordinates (1261) . i.e. x, and take values 1, 2 and 3. 

The Latin letters from the beginning of the alphabet, a, 6 , ... are used to denote the coordinates 
of 4d flat space and take values 0, 1, 2 and 3. 

The Latin letters from the middle of the alphabet, i, j, ... are used to denote the coordinates 
of the space part of 4d flat space and take values 1, 2 and 3. 

Tetrads, Second fundamental form. We partially fix the local SO(3, 1) or iS'0(4) invariance 
(which is natural to do keeping in mind the fundamental meaning of the foliation) and choose e 0 M 
to be equal to the vector normal to E^, i.e. 

e 0 M = , (28) 

where is the vector dual to the 1-form n = eNdt. Clearly, this vector is normal to the hypersurface 
t = const , and using (l27[) we see that 

/ 1 N a \ 

nn = (eAT,0,0,0) , ) and = e ' 

Then the rest of the tetrads will belong to the space tangent to Et, eT g TEj. As usual, we can 
introduce the projector on TEj 


en^ny . ( 29 ) 

The fact that h /ti , projects any vector from TM to a vector in TEt, immediately follows from that 
a ) tT = 0 and b) hffh pv = h pv . Using this, one can easily see that e,T are left invariant by h pv , 
= e^. Also, combining (l28|) . (l29|) and Zapu^v = 9pv-> one can establish that 

hp U = e ip e\ . (30) 
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Using the normal vector and the projector h /1L/ we define in the standard way a second 
fundamental form, or the extrinsic curvature, which measures how the leaves of the foliaton are 
“bended” in the ambient space-time 

K tll , = -hP\7 p n u . (31) 

From this definition and using (12911 we can establish a very useful relation: 

V fj,n v = en /J ,n rT 'V a n u - K )W . 

Covariant spin derivative. The covariant derivative on a spin bundle is defined in the usual way, 
e.g. [ 32 ]: 

= d p + uj/j, , (32) 

where = — \oo pa bl ab is a spi n connection and "f ab := ^[ 7 a , 7 6 ] are the generators of 50(3,1) or 
SO( 4). Here 7 “ are the usual flat gamma matrices, i.e. { 7 “, 7 6 } = —2rj ab , where rf b = diag(e, 1,1,1). 
The condition that covariant derivative is compatible with metric is translated into the full (i.e. with 
respect to both space-time and flat indices) covariant constancy of the tetrad: 

V p e au = d p e a i, + uj pa ebu ^' p ,u^a,p — 0 • (33) 

From here, it is easy to find the expression for oj pa b 

bJfiab — &ai/dfj,&b ^-'p.u^ap^b = p^b ; (34) 

where now is the usual space-time covariant derivative, i.e. the one acting on space-time indices 
only. 

Dirac operator. We define the standard Dirac operator without \J —1 in front: 

D = 7 ^V£ , (35) 

where 7 ^ := e a ^^ a are the curved gamma matrices, i.e. { 7 ^, 7 ^} = —2g tlLI (which is trivial by 
v = 9pu)- 

B Generalized Lichnerowicz formula 

Here we would like to calculate B 2 for a generalized Dirac operator B = 7 M V^ + F , where = 
dfi + (jO p is a metric compatible covariant derivative on a spinor bundle and F is some arbitrary 
endomorphism of this bundle. Later on we will specify the form of F. Essentially the goal is 
achieved in two steps: 1) Bring B 2 to the form B 2 = —g^d^d„ + A^d^, + B , where A and B are 
some endomorphisms of the bundle; 2 ) Using the well-known fact, see e.g. [T9] . that in this case 
there always exists a connection and an endomorphism E such that —g ,ll> d p d u + A^d p + B = 
—+ E, find these and E. 

1) B 2 = -g^d p d v + A»9 P + B 

ID> 2 = ('y fl V p + F)(YVv + F) = ^V p YVv + ^[V p ,F] + {^,F}\7 fl + F 2 = 

= -g^d p d u - (2 g^ Uv - { 7 ^, F} - - 

—g bLU {d u io p + — r a pv L0 a ) + — + 7 m (i9 m F) + { 7 ^ 0 ;^, F} + F 2 , 
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where = g au T^ v and we used the usual Lichnerowicz formula: + -j. From 

this result we identify A p and B: 

A p = -2g pv uj v + {Y,F} + ^ 

B = — g^ u (d u ujfj, + uj p uj u — r^w CT ) + — + 7 ^(d p F) + { 7 M u; M , F} + F 2 . 

2 ) B 2 = + F 

This step is standard (see e.g. nsmi). Defining 

F = F + + 

one easily establishes the desired relation: —g^d^dy + A p d p + B = — + F. In our case 

we have from (l36l) 

= 2i j -U} =: ) 

F = 7 ^(3 /i F) + { 7 ^,F} + F 2 - 

-<T K, /„} + ^'%F + F% + | . (36) 

Using that V^F = d p F + [w /t , F] and V p I u = c^/j, + [oj^YA — F^/^, the expression for F is 
simplified to become 

£ = I[7",V„F]+F 2 + I 9 „47'‘,F}{ 7 ",F} + | . (37) 

For our purposes, we need to calculate the trace of F. Clearly the trace of the commutator in 
(1371) is zero. So the only non-trivial trace is Tr(F 2 + jg p ^{ 7 ^, F}{ 7 ", F}). Using the cyclic property 
of the trace and the algebra of the 4d gamma matrices, this reduces to — Tr(F 2 — ^g /ill 'y^F r y l 'F). 
To evaluate this trace, we will take the most general form of F 

F = a + 7 p b p + 7 , (38) 

where a, b p and c pu are some functions. Because the trace of a product of an odd number of gamma 
matrices is always zero and using the well-known trace 

TKYVW) = 4 isTsT - <rg va + s pa 9 vp ) 


we get 

TV F 2 = Tr(a 2 - b^ + c^c pa 7 ^ 7 ^) = 4(a 2 - b p b p - 2 c^cT) 
g pu Tr(Y F Y F ) = g pu Tr(YYa 2 - YYYl^pbs + YY^Yl^CpaCgr,) 

= -4(4a 2 + 2b p b p ) , 

where we used [ 7 pCT , 7 ^] = —2{g av ^ p — g pu Y)- Using this one easily gets the trace of the endomor¬ 
phism F 

TrF = (-12a 2 + 8c pu c pu + R) . (39) 
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As a simple check of the result (l39l) we can easily verify that the result really should not depend 
on b fl . If F = 7 then the connection (l36l) will become + 6 ^. This is nothing but 

the connection on the twisted spinor bundle and the well-known generalization of the Lichnerowicz 
formula [33] immediately gives 

7 M V ^vt = A 6 + | + i 7 ^ , 

where Vj) = + b^, and F^y = d^b u — d v b^ is the twisting curvature. One can easily see that 

this coincides with our result (1371) for this specific form of of F. Because the trace of the commutator 
is zero, the dependence on b^ drops out of (l39l) . 

C The relation between R and R 

In this appendix we would like to obtain the relation between the Ricci scalars for a manifold 
endowed with two different Riemannian structures g and g: 

9p,u = w 

9/j.u = en^n v + ah^y . (40) 

So, we want to relate R = R[g\ and R = R\g\. Define two connections, V and V, compatible with 
g and g respectively. Then their difference is a tensor: 

~~ = —C^yUJ\ , (41) 

where uj is arbitrary 1 -form and 

C/ii/ = 2 # (^ \i9ua T — V^y) , (42) 

see, e.g. Appendix D in [34] for details. Using ([41 p in the definition of the Riemann tensor 

[V„,V„K — Rfiua^^p 

we arrive at 

V / = V / - 2 V i^ip + 2C ' 

or for the Ricci scalar 

R = rrRp, - 2 r a (y { pC\ ]a - C\^C\ ]X ) . (43) 

a) Calculating g^R^y 
From (1401) we get 

7f u = erFrF + -h^ = -g' JU + e ( 1 - - ) vFn v , 

a a \ a) 

where all the indices are raised by g^ v . Then we immediately get: 

r v Rpv = -R + e (1 - -) RpyvFrF . (44) 

a \ a 
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Because from the point of view of a foliated geometry the (3 + l)-splitting is more fundamental, we 
will re-write ()44[) in 3 + 1 form. Using the Gauss-Codazzi equations, we get: 


R^nW = K 2 - K IIU K' W + V M yn u X7 u n 11 - n^V v n v ) , 
R = ( 3 ) R+eK 2 - eK, lu K lw + 2eV, t f n A V A n^ - n^V A n A 


where K^ v 


is given by (13T1) . Combining this with (1441) we arrive at the desired result: 
r v R^v = l (3) R+eK 2 - eK^K^ + e-±^vJn x V x n» - n^V x n x ) . 


(45) 


b) Calculating - C\^ p]x ). 

From the metric compatibility of V and taking into account that g^ u = ag^ v + e(l — a)n ll n l , we 
have V p g vp = e(l — a)V M n^n p . Using this in (|42l) after some straightforward algebra we arrive at 
the following expression for C A pi ,: 


C 


A 

11V 


e(a 



— e—n u n u n p \7 p n 
a 



(46) 


From here we get 

C\^C\ ]X = . 

Because one of the indices (//, a) is always 3d and the other normal we immediately have 

rC\ [lL C\ ]x = 0 . (47) 

After some lengthy but straightforward algebra we get for the remaining term the following result: 

r a y^C\ ]a = e^-iv,KV A n A + n A V A n I/ ] . (48) 

Combining (l43l) . (l45l) . (1471) and ()48[) we arrive at the desired relation between R and R 

R = - (% +e K 2 - eK llv K llu + 2eV p f -n x ' V A n^ - n^V A n A ) . (49) 

a \a J 


D The rescaled Dirac Operator 

Here we would like to re-write the generalized Dirac operator (|5|) in the form that explicitly uses 
the re-scaled metric (1401) . i.e. we want to write D as 

D = D + S , (50) 

where D = is the Dirac operator with respect to the metric (1401) and S is some non-derivative 

part. 
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To begin with, we will re-write (j5]) using some formulas used for the (3+l)-decomposition of a 
Dirac operator (see |13| for details): 

7mV? =^D-^°K , 

+ • (51) 

Using these relations in ([5]) we get 

D = 7 °D n + ci ( % +c 2 7°/U + c 3 7 Q a a + cji + c 5 7 %“^^- = 

= (en p n v + + (y + c 2 ) 7°^ + (c 3 - ^ 7 “^+ + c 4 AT + c 5 7°7 a ^+ = 

= 7 m V^ + (y + C 2 ) 7°-^ + (c 3 - ^ 7“^- + C 4 K + C57 0 7 q ^^ , (52) 

where we used the (3+l)-decomposition of the usual Dirac operator (QJ, which immediately follows 
from (l5Tj) 

D = = 7 °D n +(3 )d_I 7 o^ + ^ 7 “^ . 

and defined the re-scaled gamma matrices: 7 ^ := ( en tJj n u + ci/i pu )'y fJ ,- Using the tetrads { e a 
adopted to the (3+1) splitting, as in the appendix [A] it easy to see that 7 ^ = ej 1 7 “ where a is 
the “flat” index and e a ^ = (n^cie^) are the re-scaled tetrads that correspond to the metric (|40D 
with a = 1/cf. The equation (1521) is not yet what we want: while the gamma matrices are already 
changed to the re-scaled ones, the covariant derivative is still defined with respect to the old metric. 
So we need to find the relation between uj p and 

Using the definition of the spin connection, = —\w l iabl ab , where u+ a ft = eai/V^e^ and the 
equation (1411) we get 

&[j,ab C a j,V fi&b /x&b "h &avC p pCj> 

a) Calculating e a v^ ^ l a l b 


e a ^V At e fe "7 a 7 b = cin^V^VV + + e^V^e/yV = 

= e ai/ V M e fe l/ 7 a 7 6 + ^<7 + - 2 ^ = 

= e av V M e b I/ 7 a 7 6 + —-— - en^SI a n v ) 7 V . 

Cl 

b) Calculating e a vC pp e b p 'y a "y b 
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Using (|46j) we get 


e au C pp e b p YY = e^-^-e ai/ (n u K pp - ec\n p n p n x \7 A n^ e b YY = 

- 2- 1 (^c\K pv e i v 'fi L - — n /J e iu n x X7xn u 'y l 'y 0 ') = 

'i V c i / 


= e- 


c 

1 -c 2 


Cl 


- (V /t! ,f+ + en p e iu n x V\n u ^ 7 V 


Combining (a) and (b) we get 


YbYY = e ay V At e b I/ 7 a 7 6 + 2-— —e^K^YY - 2e(ci - l+^+V^+V = 

Cl 

= e ay V M e b I/ 7 a 7 6 + 2 1 + 2(ci - l)n p e iv ^-YY , 

where we used the identity (B.7) from [13]: 


(53) 


n^V^+V = «+^ 7 V 


Using 


we obtain 


1 CI ~a v 0„i 1 _ d v N n i 


Y^n = ++ + —- YY K ^YY + 


= ++ + 


2 ci 

!-ci / n 




YK + Ye. pd ^- 


2 V' “ ' ' * N 

Using this result in (I52|) we finally arrive at the desired relation 


© = + ( c 2 + ^) K Y + fc 3 - 


+ c 4 iv + c 5 7 ° 7 *e i M ^^ 


(54) 


E 3+1 Metric with spherical symmetry 

The most general three-dimensional spherically symmetric metric can be written in the coordinates 
adopted to the spherical symmetry as 

^gij = diag (a 2 (r, i), r 2 6 2 (r, t),r 2 b 2 (r , t ) sin 2 (0)) . 

The spherical symmetry constrains the shift vector A* in the ADM decomposition 

ds 2 = a 2 dt 2 + ® gij(dx l + N l dt)(dx J + N^dt), 


to 

A - * = (n(r, f), 0, 0). 

Thus, the ADM decomposition of a spherically symmetrical metric can be written as 
ds 2 = (A 2 + a 2 n 2 )dt 2 + 2 a 2 ndtdr + a 2 dr 2 + r 2 b 2 dQ 2 . 
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Ill |24l[25 j it was obtained a solution to a black-hole in Aether-Einstein theory, that is also a 
solution to the low energy limit of Horava-Lifshitz gravity, 


ds 2 = f(r)dv 2 + 2 B(r)dvdr + r 2 dPl 2 , 


where v is the Eddington-Finkelstein coordinate v = t + r*, with r* such that 


dr * 
dr 


B/f, and 


fir) 

B(r) 


i + h + 

r 


Ff 


+ ... 


V 


n 

48£ r 3 
Ff r] Ff 


1 H- 1 -- H- 1 -- + 

16£ r 2 12£ r 3 


where F\ is proportional to the total mass M of a black hole as seen by a distant observer F\ = 
—2GnM and r] and £ are the parameters that enter the IR action of HL gravity (1131) . 

Switching back to Schwarzschild coordinates 


ds 2 = f{r)dt 2 + -yd 7 " 2 + r2( Fd 2 , 

we can readily identify 

N i r ) = /fir), a i r ) = -t=St > n = 0 • 

v/( r ) 

This solution implies that Earth-based experiments in our scenario would provide a contribution 
to the radial component of the Kostelecky parameter m^ in (1171) given by 


m r = 




G n M 

^2 



G 2 n M 2 t) 
4 r 2 £ 


(55) 
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